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Theoretical predictions for pp — > 2T + X cross section at -^/s = 8 TeV for the 
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INTRODUCTION 



In the recent report of the LHCb collaboration [l| a measurement of the J/ip pair pro- 
duction at pp collision energy of 8 TeV was presented. The experimental cross section value 

m; fl 

of 5.1 ± 1.1 nb appeared to be close to the predicted one in the article [2|, in which the 
leading order (LO) QCD calculations in the color singlet (CS) scheme was performed. This 
calculations take into account only the single parton scattering reactions (SPS). However, 
the double parton scattering (DPS) can also contribute to the J/ip pair production at LHC 
due to the large gluon luminosity. The cross section value of double J/ip production was 
roughly estimated within DPS mechanism is close to the SPS prediction. 

The calculation for the process gg — >■ J/ipJ/ip within LO QCD + CS + SPS is a well 
defined procedure and give a reliable results. The calculation uncertainties are determined 
by a choice of scales, parton density functions, and a mass of heavy quark. The cross section 
value within this approach can be estimated as 4 ± 2 nb Of cause there some open 
questions in this calculation: the contribution from the NLO diagrams, the influence of 
the relative quark motion in the hard part of the matrix element, and the accounting of 
transverse momenta of initial gluons. 

Recently the estimation of some high order contributions to the J/i/j pair production has 
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been done in 3]. The effects due to the transverse momenta of initial gluons was studied 

n 

in |4j within the /^-factorization approach. Also it has been shown that the quark internal 
motion could influence the cross section distribution shape 2|, as well as the cross section 

n 

value [5[. However, all these problems need further investigation. 

It is worth to note, that the standard mechanism (pQCD + CS + SPS) does not de- 
scribe well enough the invariant mass distribution of J/ip pair production lj, however it is 
show that the qualk motion in .uarkonia conld hnp.ove the agreement ff The VT and „ 
distributions for the double J/ip production were not reported yet 

In the case of double T production quarks in quarkonia have a smaller velocity 
thus the uncertainties caused by the quark internal motion are smaller than for the case of 
double J/i[) production. Therefore the formula for the cross section of the process gg — >■ 
(QQ) + (QQ) obtained within leading order of pQCD js] should describe the data on the 
double T production more precisely than the data on the double J/ip production. 

What concerns a possible DPS contribution, only simplified models can be recently used 
to estimate its contribution due to the absence of the two-particle distribution functions for 
partons in a proton. This leads to the unknown uncertainties in the cross section value. 
Using the LHCb experimental cross section value for the inclusive single J / ijj production j^] 
and the conventional value for the effective DPS cross section, the cross section value of J/ip 



pair production in the LHCb acceptance 
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llj can be estimated as 4 nb. As it was noted, 



this value is close to the result of calculation in the SPS approach. 

Later the associated production of J/ip and open charm and the double open charm 
production was measured in the same experiment [121]. It appeared that the standard mech- 
anism (pQCD + SPS) underestimates the experimental data in this channels by a factor of 
3 -T- 10. The needed cross section values can be obtained within the DPS approach, which 
implies an independent production of two cc-pairs. However it is worth to mention, that the 
DPS approach can not describe kinematical distribution shapes for the both these cases. 

In this paper we show that for the double T production in the LHCb conditions the DPS 
contribution is negligibly small, and therefore the double T production can be a better test 
for the perturbative QCD than double J/ip production. 

It worth to mention that the hadronic J/ip pair production can be researched at low 
energies in fixed target experiments, where quark- ant iquark interactions mostly contribute to 
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this process , as well as at high energies of LHC collider, where the gluon interaction 

is dominant. 

Contrary to this, the double T production can be researched only at the LHC due to the 
small cross section value (which is approximately 10 3 times smaller than for the J/ip pair 
production), as well as the large invariant mass needed to produce an T pair. It could be 
expected, that the gluon fusion gives the main contribution into the process at the LHC. 

In this work we are mainly interested in the low invariant mass region as it is most 
accessible in the LHCb experiment, which specializes in heavy quarks studies. We assume, 
that the next leading order corrections are small for the LHCb experiment conditions. It 
should be noted here that at large invariant masses the high order corrections change the 
cross section behaviour and it becomes a constant, as well for the double J/^-production [161 ]. 
while within the LO of perturbative QCD the cross section for these processes decreases with 
gluonic energy as ~ 1/s 2 , as it is shown in [8j. 

ATLAS is another LHC experiment, where the T pair production can be researched. 
Despite the pt cuts applied for muons from T decays, ATLAS can reconstruct T-mesons 
without cuts on their pt jl3] due to the large energy release in the leptonic T decay. The 
efficiency and the total number of reconstructed events in the low-p^ region is not expected 
to exceed the LHCb values, however the cross section of double T production in the ATLAS 
experiment is predicted to be larger than in LHCb. 

It is worth to note that in the leading order of a s the quarkonia pair production obeys 
selection rules, which are similar to those for quarkonia decays. Selection rules in the 
charmonia decays are confirmed by measurements of J/ip and \ c meson widths. In the 
production case two initial gluons in a color-singlet state are C-even. This is why the 
production of T-, r\y- or xv meson pairs is allowed, while the combined production of two 
particles having a different C-parity (such as T + rj b and T + Xb) is prohibited. 

These selection rules put restrictions on the feed-down from the higher states. For instance 
the feed-down from the T + Xb production is expected to be less than from the Xb + Xb one, 
because the production of the former final state is forbidden in the leading order of pQCD. 
Meanwhile the feed-down from the T + T' channel can be quite large. 



1 the simultaneous production of two J/T/j-mesons was first observed in f982 by the NA3 collaboration in the multi- 
muon events in pion-platinum interactions at f 50 and 280 GeV and later at 400 GeV in proton-platinum collisions. 
These data is fairly described within the CS model 
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We expect that the color-octet mechanism negligibly contributes to the double T produc- 
tion. Indeed, the color-octet contribution to the T production should be essentially smaller 
than the color-octet contribution to the J / ip production, because the color-octet matrix el- 
ement is suppressed by the second power of relative velocity of the quarks inside the meson. 
But even for the doubly J/tp production at LHCb the color-octet should not be taken into 
account 2 . 

The T-meson pair production will be investigated in the decay mode T — > fi + fi~ for 
the both quarkonia: 2T — > fi + fi~ fi + fi~ . The region of invariant masses of four muons 
near the threshold of the T-pair production can be interesting due to the opportunity for 
two diquarks [6% + \bb\ to form a bound state — a tetraquark, which could decay to 
T + — > fi + [i + fi~ . An attraction between 3 and 3 color states does not exclude 
such a possibility, especially since similar exotic states like the Y~(3940) resonance, decaying 
to J/ipu {j^J, and the X(4140) resonance, decaying to J/ijjcf) 21, 22], have recently been 
observed. It should be mentioned that according to the LHCb experimental data, there is 
no evidence of the 4c-tetraquark, decaying into two J/ip mesons. However [66] g is essentially 
smaller than [cc]g and therefore the model of two attracting doubly heavy diquark can be 
more applicable for the 46-tetraquark. 

The second section of our article is devoted to the non-resonant production of T-meson 
pairs in the gluon-gluon interaction. In the third section the cross section of this process 
at the LHC energy of 8 TeV is calculated taking different experimental restrictions into 
account. The fourth section the production of 46-tetraquark at ther LHCb experiment is 
discussed. 

II. DOUBLE T PRODUCTION 

The standard method to describe the double T production is based on following assump- 
tions: 

1. the LO QCD can be used to calculate a matrix element of four heavy quark production; 

2. b and b quarks in T are in a color singlet state; 



2 As it is shown in [la . Ila ] the octet contribution becomes significant only at transverse momentum pr > 15 GeV, 
which corresponds to a big invariant mass of the J/tp pair. 
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3. the internal motion of quark in the meson do not essentially influence the cross section 
value and the distribution shapes (^-approximation); 

4. the both T mesons are produced in a single gluon-gluon interaction (a single parton 
scattering, SPS). 

In the leading order of perturbative QCD there are 31 Feynman diagrams describing a 
color-singlet T pairs production in the gluonic interaction. 

We do not consider the contribution of quark- ant iquark interaction, which is negligibly 
small at LHC energies. 

The hadronization of 66-pair into the final T-meson within the ^-approximation is ac- 
counted for by the wave function of this particle at origin: 

V'T(is)(r)| r=0 = 0.635 GeV 3 / 2 , 

V>Wr)| r=0 = 0.455 GeV 3 / 2 , (1) 
V>T(3S)(r)| r=0 = 0.400 GeV 3/2 . 

These are only nonperturbative parameters in the matrix element of gg — > 2T process. The 
values (CQ) are extracted from the leptonic widths of T-mesons, neglecting QCD corrections, 
as we do not take these corrections into account in our matrix element. 

As every calculation, our calculation of T pair production is affected by the scale 
uncertainty. This one is of the same order of magnitude, as the 6-quark mass uncertainty. 
Unlike the scale choice, the choice of the 6-quark mass also affects on a relative yield of 
different bottomonia pairs. 

The conventional b mass value in the matrix element equals a half mass of produced 
bottomonium. Alternatively, one can choose one and the same nib value for all bottomonia 
pairs. We present the results obtained for the both these cases. For the second one we use 
nib — 4.73 GeV in the matrix element while the mass of the particular T-meson is taken in 
the phase space. 

To estimate the feed-down from higher excitations we use the PDG |23] branching frac- 
tions for T(2S) — > T(1S)+X and T(3S) — > T(1S, 2S)+X decays. The impact of feed-down 
on the cross section values is less than the uncertainties owing to the scale and the nib choices. 
Nevertheless, it should be taken into account, because it effects on the relative yields of dif- 
ferent bottomonia pairs. The T(1S')T(1S') pair production is most influenced by feed-down, 
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Mode 


m& = tot/2 


m b = 4.73 GeV 


without feed-down 


with feed-down 


without feed-down 


with feed-down 


T(1S)T(1S) 


36 pb 


45 pb 


36 pb 


48 pb 


T(2S)T(2S*) 


5.3 pb 


6.0 pb 


8.6 pb 


9.8 pb 


T(3S)T(35) 


2.3 pb 


2.3 pb 


4.9 pb 


4.8 pb 


T(15)T(25) 


27 pb 


33 pb 


35 pb 


44 pb 


T(15)T(35) 


18 pb 


20 pb 


27 pb 


31 pb 


T(2S)T(35) 


7.0 pb 


7.4 pb 


13 pb 


14 pb 



TABLE I: The double bottomonium production cross sections in 8 TeV pp-collisions. No kinematical cuts 
are imposed. A relative uncertainty of 1.4 due to the scale choice is assumed. 

where it leads to the 30% increase of the yield. 

The transverse momenta of the initial gluons are taken into account within the Pythia 
machinery. 

In Tab. [J we present the double bottomonium production cross sections (without kine- 



It 



rese predictions were done 



24j at the scale equal to an 



matical cuts) estimated for the different parameters values. A 
for the 8 TeV pp collisions using the CTEQ6LL p.d.f. set 
averaged transverse mass of quarkonia. The a s in the matrix elements of hard subprocess 
is also taken at this scale. 

The corresponding cross section distributions on the invariant mass are shown in Fig. [TJ 
The distribution shapes do not depend essentially on the hard scale choice, while the cross 
section value changes by factor of approximately 1.4 with the variation of scale by a factor 
of 2. 

Despite the difference in masses and in feed-down contributions, shapes of distributions 
over the transverse momenta and rapidity of different T mesons do not differ significantly. 
At least the uncertainties caused by ambiguity in scale selection are larger. This is why, in 
what follows, we discuss kinematic properties of T(IS') pair production only. 

Here we consider the T pair production at LHCb and ATLAS experiments. 

The LHCb experiment, being designed for the 6-physics studies, allows to measure the 
T pair production without px cutoff. Contrary to this, as a rule, the ATLAS experiment 
allows to study only relatively large px- However, namely the T production is an exclusion. 
Contrary to the J/ip meson leptonic decay, the energy release in the T decay is large. This 
allows ATLAS to measure T mesons even at low pr- The production cross section in ATLAS 
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FIG. 1: The distributions over the invariant mass of T-meson pairs. A half of the meson mass is used for the 
nib in the left plot and fixed value of mj, = 4.73 GeV is used in the right. From left to right curves correspond 
to the T(15)T(15), T(1S)T(2S), T(1S)T(3S), T(25)T(2S*), T(25)T(3S) and T(3S*)T(35) production. 
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FIG. 2: The rapidity distribution of the T-pair in the T-pair production process with the LHCb rapidity 
cut. 



conditions appears to be even larger than in the LHCb due to the wider rapidity window 
(-2.5 < y < 2.5 instead of 2 < y < 4.5 for LHCb). 

Within the SPS approach we predict the hump in the rapidity distribution of single T 
meson in the double T production at LHCb, as well as the hump in the rapidity distribution 
of single J / if) meson in the double J / if) production. This hump is caused by the LHCb cutoff 
on rapidity, which leads to the strong rapidity correlation (see Fig. [2]). 
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FIG. 3: Distribution over the rapidity of single T in the T-pair production process at ATLAS. 



Mode 


ATLAS 


LHCb 


mi, = my /2 


m b = 4.73 GeV 


mi, = rot/2 


m b = 4.73 GeV 


T(1S)T(1£) 


31 pb 


33 pb 


6.2 pb 


6.6 pb 


T(2S')T(25) 


4.2 pb 


6.8 pb 


0.8 pb 


1.3 pb 


T(35')T(35) 


1.6 pb 


3.4 pb 


0.3 pb 


0.6 pb 


T(1S)T(2£) 


23 pb 


30 pb 


4.5 pb 


6.0 pb 


T(1S)T(35) 


14 pb 


21 pb 


2.8 pb 


4.2 pb 


T(2S)T(35) 


5.2 pb 


9.6 pb 


1.0 pb 


1.9 pb 



TABLE II: The cross sections of doubly bottomonium production in 8 TeV pp-collisions. ATLAS and LHCb 
kinematical cuts are imposed. 

The rapidity distribution in ATLAS conditions is presented in Fig. |3j 
Despite the essential difference in the kinematical conditions at these facilities the dis- 
tributions over the pt of single T in the pair production and over the pt of T pair have 
practically the same shapes. This why we perform here these distributions without rapid- 
ity cuts (see Fig. H] and Fig. [5]). As one can see, the different choices of hard scale value 
significantly affect the slope of both distributions in the high pt area. 
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FIG. 4: The pr distributions of the single T in the T(15)T(15) production process at LHC depending on 
the scale choice. 
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FIG. 5: Distribution over the pr of the T-pair in the T(15)T(15) production process at LHC depending 
on the scale choice. 



In the LHC conditions a huge density of \ow-x gluons leads to the increase of multiple 
gluon-gluon interactions probability within a single proton-proton collision. In the DPS 
approach, which implies a particles production in two independent subprocesses, the cross 
section is written down as follows: 



A B 
AB _ m °SPS°SPS 
°bPS — ~7> ~ • 
^ Ocff 



(2) 
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where the parameter a c g = 14.5 mb was measured in the four jets and three jets plus photon 
modes by the CDF and DO detectors Q, Q, and the parameter m equals 1 for identical 
subprocesses and 2 for different ones. As it was shown in 25j-|27j, the DPS mechanism could 



play an essential role in the double J/ip production, because the predictions for the double 
J/ip production within SPS and DPS have the same order of magnitude. For the double T 
production the DPS mechanism contributes only about 10% of the SPS cross section and, 
therefore, can be neglected (see [271]) . 



III. THE POSSIBILITY OF OBSERVATION OF THE 46-TETRAQUARK AT LHCB 



Within QCD two heavy quarks in the 3 color state attracts each other, forming the 
[QQ]s diquark. An attraction between 3 and 3 color states does not exclude a possibility 
to observe the heavy tetraquark [QQMQQJs- The spectroscopy of such systems can be 
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. Within 



easily investigated under assumption that these diquarks are almost point-like 
this approach it was predicted that all 46-tetraquark states have masses below the T + T 
threshold. Therefore such states can not contribute to the double T production. Nevertheless 
the 46 tetraquark could be searched in the decay mode [56] § [66] 3 — > T + T* — > Tji + n~ . 
According to 2^] the ground state of 46 tetraquark splits into tree states due to the spin- 



spin interaction: 



++ 
1+- 
2++ 



M = 18.754 GeV, 
M = 18.808 GeV, 
M = 18.916 GeV, 



M-M th 
M-M th 
M-Mth 



-544. MeV, 
-490. MeV, 
-382. MeV. 



It worth to note, that at the present time there are no reliable methods to estimate the 
cross section value of the 46 tetraquark production. 

IV. CONCLUSIONS 

The observation of J/ip pair production has opened an interesting discussion about the 
contributions of SPS and DPS mechanisms into the production process. The SPS contri- 
bution is estimated in the framework of well known and well tested pQCD approach. The 
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second mechanism, DPS, is not well studied experimentally yet. However it allows to under- 
stand the large cross sections of the J/ip plus open charm production and the double open 
charm production at LHCb. For the double J/ip production the SPS and DPS approaches 
give comparable predictions, which are of the same order with the experimental results. 

The situation is different for the double T production. The both SPS and DPS calcu- 
lations for it stand on the same basis as for the double J/ip production, but for this case 
the SPS prediction is about an order of magnitude larger than the DPS one. Thus, the 
study of T pair production will allow us to understand, if we describe the SPS production 
of quarkonia correctly. 

The observation of double J/ip- meson production [lj was reported by LHCb soon after 
the measurement of inclusive J/i[) production |9j]. The measurement of T-meson inclusive 



production have been recently reported in 



291 ] . It is based on the 25 pb 1 integrated lu- 



minosity collected in 2010. The integrated luminosity collected both in 2011 and in 2012 
exceeds 1 fb _1 . This corresponds to more than 40 times lager statistics. Taking into account 
the estimated here cross section value and the leptonic branching of T-meson (2.5%) one 
can expect several dozens of events with two T mesons in both 2011 and 2012 data. 
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